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Abstract 

We report on a search for right-handed W bosons (Wi?). We used data col- 
lected with the D0 detector at the Fermilab Tevatron pp collider at y/s =1.8 
TeV to search for Wr decays into an electron and a massive right-handed 
neutrino — > e^Npt- Using the inclusive electron data, we set mass limits 
independent of the Nr decay: mpi/^ > 650 GeV/c^ and mwj^ > 720 GeV/c^ 
at the 95% confidence level, valid for mNj^ < ^my/^ and m^^ <^ rnwa 
spectively. The latter also represents a new lower limit on the mass of a heavy 
left-handed W boson (W') decaying into ez/. In addition, limits on niWjt '^^lid 
for larger values of the mass are obtained assuming that Nr decays to an 
electron and two jets. 

Right-handed W gauge bosons (Wr) are additional intermediate vector particles that 
arise in extensions of the Standard Model (SM) such as the left-right symmetric model 
(LRM) [|I|. In the LRM, an enlarged SU{2)r x SU{2)l x f/(l) symmetry group replaces 
the SU{2)l X U{1) group of the SM. As a result of the additional symmetry, three new 
gauge bosons, two charged and one neutral Z', appear along with massive right-handed 
neutrinos (Nr). 

5 



In this Letter, a direct search for Wr bosons with mass greater than 200 GeV/c^ which 
decay into an electron (or positron) and a massive right-handed neutrino, Wn — > eNji 
is reported. The Nfj, is assumed to decay promptly through the right-handed charged 
current into a mode that depends on the mixing angle ^ between Wl and Wr. If the 
mixing is negligible (no mixing case), the Nr will decay into an electron and an off-shell 
Wr, Nr — > eW^. The right-handed neutrinos from other lepton families are assumed to 
be at least as massive as the electro n-A?"/?. Therefore, the off-shell Wr can decay only into 
quarks, Wr qq' . On the other hand, if the mixing is large, the Nr decays into an electron 
and a W boson, which decays into quarks two thirds of the time. In both cases the decay 
chain leads predominantly to a final state with eeqq. 

Previous direct searches at hadron colliders yielded the lower limits > 261 GeV/c^ 
0, valid for any value of the mass of the right-handed neutrino, and niwji > 652 GeV/c^ [||, 
valid only for a light right-handed neutrino (mTv^ -C mwji) that does not decay or interact 
within the detector. Indirect searches based on low energy phenomena such as fi decay, 
the Kl-Ks mass difference, and neutrinoless double beta decay provide additional stringent 
lower limits [Q. Limits from direct and indirect searches depend, however, on the assumed 
values of the elements of the mixing matrix for the right-handed quarks, the coupling 
constant Qr, the mass and type (Dirac or Majorana) of the right-handed neutrinos, and the 
mixing angle ^. The most general limit is ■ > 300 GeV/c^ ||^. 

Two different methods, corresponding to different values of the ratio Rm = f^Nn/^WR, 
are used for this search. For < |, the products of the Nr decay are not likely to be 
well separated, making their individual identification difficult. Therefore, the transverse 
momentum spectrum of the Wr decay electron, which is expected to be hard and to have 
a distinctive Jacobian peak at (m^^ — mj^^ /2mw^, is used as a signature. A search for 
such a peak, henceforth referred to as the peak search, is carried out using the high-py 
inclusive electron data. This method does not discriminate between helicities of the W boson. 
Therefore, the peak search is also sensitive to heavy left-handed W bosons {W) which decay 
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into an electron and an electron neutrino W ev. For Rm > |, the products of the Nr 
decay are likely to be well separated, making possible the detection of the exclusive final state 
with two electrons and two jets. After requiring the two electrons to be inconsistent with 
Z ^ ee decay, the background due to other known physics processes is small. Therefore, a 
simple counting experiment, referred to here as the eejj search, is performed. The analysis 
presented here is based on approximately 79 pb^^ of data collected during two Fermilab 
Tevatron pp collider runs at a/s =1.8 TeV from August 1992 to May 1995. 

The D0 detector consists of three major subsystems: a central tracking system with no 
magnetic field, a hermetic uranium-liquid argon sampling calorimeter, and a muon magnetic 
spectrometer. The calorimeter has fine longitudinal and transverse segmentation in pseudo- 
rapidity {rf) and azimuth (0) that allows electromagnetic showers to be distinguished from 
jets. It provides full coverage for |?7| < 4 with energy resolution 15%/y^-E'(GeV) for elec- 
tromagnetic showers and 80%/ \J E{GeY) for hadronic jets. The central and forward drift 
chambers are used to identify charged tracks for \ri\ < 3.1 and to locate the primary vertex. 
A more detailed description of the D0 detector can be found elsewhere 

To identify electrons 0], the presence of an isolated electromagnetic energy cluster with 
shape consistent with that of an electron (as determined from test beam measurements) is 
required. In addition, an associated charged track that matches the calorimeter cluster in 
7] and (f) and with an ionization in the drift chambers consistent with that of a minimum 
ionizing particle must be found. Jets are reconstructed using a cone algorithm with a cone 
radius of 0.5 in ?7-0 space. 

For the peak search, events were collected using a single electromagnetic cluster trigger. 
Offline, the inclusive high-p^ electron events were selected by requiring an electron candidate 
with > 55 GeV/c and |?7e| < l-l- To reduce the multijet background (QCD) from events 
with a jet misidentified as an electron, strict electron identification criteria were imposed. 
The 101 events with > 100 GeV/c were scanned to search for anomalies; we discarded 
one event which was consistent with being a high energy cosmic ray muon that showered in 
the electromagnetic part of the calorimeter, mimicking an electron. 
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FIG. 1. (a) Electron transverse momentum and (b) transverse mass (formed by the electron 
and ^2") distributions of the inclusive high-py electron sample. 

The primary background in the peak search is due to highly off-shell and large-p^ W and 
Z boson production. These processes were simulated using a Monte Carlo (MC) program 
based on a theoretical calculation of the bosons' pr [§ and on the bosons' line shape obtained 
using the PYTHIA MC program, with a simple detector simulation. The QCD background 
was modeled using the collider data. 

A simultaneous fit to the transverse mass (m-r) distribution, formed by the electron and 
the missing transverse energy and to the electron transverse momentum (p^) distribution 
was performed. A binned maximum likelihood fit was used to find the contributions of the 
combined W and Z boson backgrounds and the QCD background [jlOl . Figure ^ shows the 
p^ and rriT distributions with their corresponding fits. The confidence level (CL) is 71% for 
the p^ fit and 90% for the fit. The presence of Wr cNji decays would appear as an 
excess in a few consecutive bins in the p^ distribution. No evidence for such an excess is 
observed. 

The acceptance and p^ distribution of the signal were obtained for a grid of points in 
the (mvy^,mAr^) plane using PYTHIA MC samples with a detector simulation based on the 
GEANT program . The 95% CL upper limit on the number of Wr events was obtained 
by integrating the probability of the presence of a Wr component in the measured p^ 
distribution for every point in the grid. This was converted into an upper limit on the 
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FIG. 2. 95% CL upper limit on aB as a function of the Wpt boson mass. Limits are shown for 
three values of the A'^^ mass. 

cross section times branching fraction (crB) by normalizing to the measured W and Z boson 
production cross sections [|1^] using the observed W/ Z component in the initial simultaneous 

and rriT fit and the acceptances as calculated from MC simulation. 

The resulting background subtracted upper limit, including the effect of systematic un- 
certainties (dominated by a 7.6% uncertainty in the W/Z background normalization), is 
shown in Fig. Also shown is a second order theoretical calculation []T^ of cxB assum- 
ing Qji = ql and = V^. The next to leading order MRS(H) [|14| parton distributions were 
used for the calculation. The branching fraction B{Wr —>■ cNr) was calculated taking into 
account the Nr and t-quark masses and assuming ■mN<^ = rriiy^ = rriN-^- For small mass, 
this fraction approaches the naive ^ value. Figure ^ shows the corresponding excluded mass 



region. The contours are shown for different values of the LRM parameters gp; and [|T5 
The extreme effect of varying is illustrated by displaying the contour for a mixing matrix 
with = 1 (thus Vj^ = for unitary), suppressing the primary ud — > Wji production 
mechanism. Because the limit from this part of the search was extracted from the inclusive 
distribution, without additional topological requirements, it is valid irrespective of the 
specific decay mechanism for the Nji or the W helicity. 
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FIG. 3. 95% CL excluded Wr mass region from the peak search. The Unes represent the 
contours for different values of the LRM parameters. The diagonal line is the kinematic limit for 
the Wr eNfi decay. 

For the eejj search, events were selected using a trigger that required two electromagnetic 
energy clusters, each with Ej- > 20 GeV. After event reconstruction, 22 events had two 
good isolated electrons with Et > 25 GeV and two or more jets with Et > 25 GeV 
within a pseudorapidity range \r]e,j\ < 2.5. Events consistent with Z+jets production were 
rejected by demanding that the invariant mass of the two electrons rriee be outside the range 
70 < rriee < 110 GeV/c^. Two events remained in the sample and were therefore considered 
Wr candidates. 

The largest background to the eejj signal is multijet production (QCD) with two jets 
misidentified as electrons. To calculate this background, the invariant mass spectrum of 
the jet pair with the largest electromagnetic fraction in events with four or more jets was 
found. This distribution was then scaled by a factor determined from a two-component fit 
to the inclusive dielectron data using the dielectron invariant mass spectrum from Z, 7* MC 
and the measured inclusive dijet invariant mass spectrum. The background from Z, 7*+jets 
production was estimated by scaling the number of observed events in the peak of the rriee 
distribution, in events with two or more additional jets, by the tail-to-peak ratios obtained 
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TABLE I. Background estimates and event yields for the eej and eejj samples. 



Background 


Event Yield for 79.0 ± 4.3 pb" 


-1 


Process 


eej 




eejj 




12.84 ± 2.31 




1.26 ± 0.34 


tt 


0.61 ± 0.35 




0.43 ± 0.16 


WW 


0.13 ± 0.02 




0.01 ± 0.01 


QCD 


9.90 ± 4.01 




1.38 ± 0.68 


Total 


23.48 ± 4.64 


3.08 ± 0.78 


Observed 


23 


2 



from MC. Additional background is due to tt and WW production. The yield from tt was 
obtained using a Monte Carlo sample with a detailed detector simulation and the measured 
6.4 ± 2.2 pb [|I^ cross section. For the WW background, a sample of MC events and the 
theoretical cross section were used. To verify the background estimation, the yield of the 
above processes to a final state with two electrons and one or more jets was also calculated. 
The background estimates and event yields are summarized in Table | for the eej and eejj 
final states. 

As for the peak search, the signal acceptance for the eejj search was calculated for a 
grid of points in the {rriw^^m]^^) plane using MC simulation. The electron identification 
efficiency was determined from Z — > ee data. Example signal efficiencies for the no mixing 
case are (15.0 ± 1.7), (10.1 ± 1.4) and (1.0 ± 0.4)% for {mw^^niN^) = (650,200), (400,350) 
and (400,50) GeV/c^, respectively. For the large mixing case the corresponding efficiencies 
are lower due to the smaller Nji — > eqq branching fraction. Also, for the large mixing case 
the search was restricted to rrif^^ > 90 GeV/c^ since the efficiencies vanish when ttinj^ ~ mw 
due to a threshold effect. 

Given no observed excess of events beyond the expected background, we set a 95% CL 
upper limit on aB using a Bayesian approach [jl^ with a flat prior distribution for the signal 
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FIG. 4. 95% CL excluded region of Wr mass from the eejj search for the no mixing case. 

cross section. The uncertainties on the overall efficiency (10-20%), the integrated luminosity 
(5.5%), and the background estimation (25%) were included in the limit calculation with 
Gaussian prior distributions. The resulting background subtracted upper limit is plotted in 
Fig. 1^, while Fig. H shows the excluded region of the (mvy^,mAr^) plane for the no mixing 
case. 

In conclusion, no evidence for the production of right-handed W bosons was found. From 
a peak search we set mass limits independent of the Nji decay: mvi/^ > 650 GeV/c^ and 
^Wr > 720 GeV/c^ at the 95% CL, valid for mjy^ < \mwj^ and mMp. -C respectively, 
assuming SM coupling {qr = Ql and = V^). Also from the peak search, we set a mass 
limit of miyi > 720 GeV/c^ at the 95% CL, extending the previous most stringent limit 
for heavy left-handed W bosons which decay into eu. In addition, limits on valid 
for larger values of the mass were obtained assuming that the Nji decays to an electron 
and two jets. Figure ^ summarizes the results of the two methods used for the search as 
an exclusion region in the {mwR,^^^) plane. These limits on place stringent, though 
model dependent, limits on possible V + A couplings. 

We thank the Fermilab Accelerator, Computing, and Research Divisions, and the support 
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